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Abstract 

We derive the theory of the quantum (zero temperature) superconductor to metal transition in 
disordered materials when the resistance of the normal metal near criticality is small compared to 
the quantum of resistivity. This can occur most readily in situations in which "Anderson's theorem" 
does not apply. We explicitly study the transition in superconductor-metal composites, in an s- 
wave superconducting film in the presence of a magnetic field, and in a low temperature disordered 
d-wave superconductor. Near the point of the transition, the distribution of the superconducting 
order parameter is highly inhomogeneous. To describe this situation we employ a procedure which 
is similar to that introduced by Mott for description of the temperature dependence of the variable 
range hopping conduction. As the system approaches the point of the transition from the metal 
to the superconductor, the conductivity of the system diverges, and the Wiedemann-Franz law is 
violated. In the case of d-wave (or other exotic) superconductors we predict the existence of (at 
least) two sequential transitions as a function of increasing disorder: a d-wave to s-wave, and then 
an s-wave to metal transition. 



1 



I. INTRODUCTION 



The quantum (zero temperature) transition from a superconducting to a non- 
superconducting ground state is the poster child of quantum phase transitions. This tran- 
sition is induced by changing external parameters at zero temperature T. 

In this article we consider three representative problems in which a direct quantum phase 
transition occurs from a superconducting to a metallic phase in which kpl ^> 1: the case of 
a composite of a superconducting and a non-superconducting metal in which the effective 
interaction between electrons changes sign as a function of position (Section II), the case 
when s-wave superconductivity is destroyed by an external magnetic field (Section III), and 
the case when d-wave (or other exotic) superconductivity is destroyed by quenched disorder 
(Section IV). Here kp is the Fermi wave- vector, and I is the electron mean free path. 

In three dimensions (3D), there is no question concerning the existence of both a super- 
conducting and a metallic phase. In 2D, the existence of a metallic phase is problematic. 
However, for kpl ^> 1, single particle localization occurs on such large length scales that 
its effects are mostly unobservable. Therefore, for most of this article, we will ignore the 
fundamental, but for our purposes purely academic question of whether or not a 2D inter- 
acting system of electrons can ever exhibit a metallic state in the asymptotic limit of zero 
temperature and infinite volume. 

Before proceeding to discuss the findings of the present study, it is worth commenting 
briefly on the broader context. The central insight underlying the modern theory of critical 
phenomena is that, due to the divergent correlation length at criticality, the properties of 
the critical state are "universal." An important extension of this is the idea that, in sys- 
tems with quenched disorder, the variations in local environments are self-averaging, so a 
near-critical system can be treated in terms of an effective, translationally invariant field 
theory While this approach has had notable successes for the theory of classical 

phase transitions, it is more problematic in the case of quantum critical phenomena. This 

n n n 

is most dramatically illustrated by the case[3J, U, |5J of the quantum critical point in the 
random transverse-field Ising model, where the physics of "rare events" results in the exis- 
tence of a "quantum Griffith phase" in which, for a finite interval of parameters including 
the critical point, the susceptibility diverges as T — * 0. For somewhat analogous reasons, 
a generic feature that characterizes the transitions in all three cases mentioned above is 
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that, at criticality, the spatial distribution of the superconducting order parameter is highly 
inhomogeneous. It is concentrated in "superconducting puddles" where, due to randomness, 
superconducting order is locally anomalously favorable, and the distance between "opti- 
mal" puddles is parametrically large. The transition occurs when the Josephson coupling 
between optimal puddles (which falls with a power of the separation) times the exponen- 
tially large local superconducting susceptibility on a puddle is strong enough to stabilize a 
macroscopically phase coherent state. 

Near enough to the quantum phase transition and at low enough temperatures, where 
the correlation length is large compared to the distance between superconducting puddles, 
there is still, presumably, universal behavior described by appropriate critical exponents. 
However, a consequence of the large distance between optimal puddles is that this universal 
quantum critical regime is parametrically narrow. Conversely, there is an anomalously broad 
portion of the phase diagram in which the correlation length is comparable to or smaller than 
the distance between optimal puddles, but large compared to the size of an optimal puddle, 
where quantum fluctuations of the superconducting order-parameter dominate much of the 
physics. This broad quantum but not quantum critical regime is one of the characteristic 
signatures of the inhomogeneous nature of the critical state. 

It is obvious that the presence of significant superconducting correlations in the "metallic" 
state near to the superconductor to metal critical point makes it highly anomalous: Its zero 
temperature conductivity diverges at the point of the transition, and can be much larger 
than the Drude value "near" criticality. In the same regime, the Hall conductivity decreases 
with respect to the Drude value and vanishes at the point of the transition. The Wiedemann- 
Franz law in such metals is also clearly violated. 

One remarkable implication of the present analysis is that, in the case when d-wave 
superconductivity is destroyed by disorder, there are (at least) two quantum transitions: 
the first from a globally d-wave to a globally s-wave (although, possibly, still locally d-wave) 
state, and the second to the "normal" metal. Another outcome of this picture is peculiar 
temperature dependencies of the physical parameters of the near critical superconductor. 
Some of our findings are summarized in the schematic phase diagrams shown in Figs. 1-3. 

To conclude this introduction, we would like to discuss the relation between our paper 

and those BBS 

in which it has been proposed that the quantum transition, especially 
in 2D, takes place between the superconducting and an insulating state. Several lines of 
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reasoning led to the inference that near the T = quantum critical point fcjri ~ 1 , in which 
case localization (which we neglect) would necessarily be a serious issue: 

1) Where Anderson's theorem applies at the level of mean-field theory, such as in the 
case in which s-wave superconductivity is destroyed by increasing disorder, the localization 
length must be comparable to or shorter than the coherence length in order for the disorder 
to have any substantial effect on T c , at all 9|. 

2) It has been shown that in a system of superconducting grains linked by resist iyely 
shunted Josephson junctions, quantum fluctuations of the order parameter llil, UM 
are strongly suppressed so long as G^^ ^> 1, independently of the strength of the Josephson 
coupling between puddles! Here G^ e ^' is a dimensionless shunt conductance measured in 
units e 2 /h. An apparent implication of this result is that, so long as a small portion of a 
highly conducting system is superconducting, at low enough temperature the system will 
achieve global phase coherence so long as the dimensionless effective conductance, G^ e ^\ is 
large compared to 1. 
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of 2D bosonic models of the transition 



3) It was found in several theoretical studies 
that there is a universal value of G^^^ = G c ~ 0{1) at the point of superconductor- 
metal transition. (These models assume the absence of gapless quasiparticle excitations, 
and therefore ignore dissipation of the sort that is represented by the shunt resistance in the 
previously discussed models.) 

4) A large portion of the experimental realizations of such transitions involve two di- 
mensional (2D) systems, such as films. Here, the transitions are often referred to as a 



"superconductor to insulator transition" [6| on the basis of the widely held theoretical belief 
that metallic phases are forbidden in 2D due to single particle localization [k], [14J] - any 
non-superconducting phase is thus expected be insulating at zero-temperature. 

Concerning point 1), in the present paper, for the most part, we consider problems in 
which Anderson's theorem does not apply, either due to the symmetry of the order parameter 
{e.g. d-wave) or the breaking of time reversal symmetry {e.g. by an applied magnetic field). 
Concerning point 3), we consider cases in which gapless quasiparticles are present near 
criticality, so the applicability of a bosonic model is questionable. Concerning point 4), as 
mentioned above, 2D localization is negligible, and hence a non-superconducting phase is 
"effectively metallic," whenever the parameter kpl is sufficiently large, since in this limit, 
the localization length is exponentially large, £; oc ~ Zexp(|G2L>)- Here G^d oc k 2 F ld ^> 1 is 
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the dimensionless conductance per square measured in units of e 2 /h of a 2D film of thickness 
d. (For a review, see 151].) Finally, concerning point 2), a large portion of the discussion 
in the present paper follows from the same considerations. The differences between the 
present results and those of the earlier studies spring from the fact that the effective model 
we develop from microscopies differs in a subtle manner from the phenomenological Ohmic 
heat bath considered in those earlier studies. This difference permits a transition to a phase 
incoherent state even under conditions in which G^f^ 1; however, a residual consequence 
of the same physics discovered in those earlier studies is that this transition occurs when 
the superconducting puddles are extremely dilute and so are weakly Josephson coupled to 
one another. This is one of the central results of the present study. 



A. Effective action for the quantum superconductor-metal transition 

In the present section, we develop the general features of the effective actions that govern 
the quantum fluctuations of the order parameter near criticality. Formally, such effective 
actions are obtained by integrating out the fermionic degrees of freedom, and all high energy 
collective modes, leaving us with a set of degrees of freedom, Aj = |Aj| exp[i0j], identified 
as the phase and modulus of the order parameter on puddle i. In detail, the various terms 
are sensitive to the specific physical circumstances, but the overall structure of the effective 
action is the same in all cases studied in the present paper. 

In the case of small puddles embedded in a normal metal, where the value of the order 
parameter is small, the Andreev reflection of quasiparticles from the metal-superconductor 
boundary is ineffective and can be taken into account in perturbation theory. Then, the 
imaginary time effective action that governs the superconducting fluctuations near criticality 
is of the form 

S = Y1 | a i/ 

+ |/ ^ ( g:y } + /^}] + - 

Hj[{A}] = -(1/2) X^A'Aj ■+ c.c] (1) 

where r is imaginary time, j labels the randomly distributed superconducting puddles, 7 
is the parameter that tunes the phase transition (e.g. the magnetic field in units of the 
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critical magnetic field), jj c is the critical value of 7 in the j t h puddle , A is the magnitude 
of the order parameter deep in the superconducting state , aj and (3j depend on the local 
structure of the superconducting puddle (as discussed below), and is the Josephson 
coupling between nearby puddles. The . . . symbol represents high order terms in A that 
are negligible at the phase transition, including non-local quartic terms involving the order 
parameter on more than one puddle. (Some representative aspects of the derivation of Eq. [I] 
are sketched in the Appendix.) 

The first and the third terms in Eq. [1] reflect the dynamics of the BCS Cooper instability 
3,0, 3, and hence 

Oi^uVi, A~^/A (2) 

where Vi is the volume of the i th puddle and v is the density of states in the surrounding 
metal. (In two dimensional cases, naturally, Vi is the area of the puddle and v = u^jjd is the 
two dimensional density of states.) To illustrate this, set Jij = and consider the dynamics 
of small amplitude fluctuations of the order parameter on an isolated puddle that is on the 
verge of becoming superconducting (1 ^> (7^ — 7) > 0) as described by the linearized version 
of Eq. Q] 

Comparing this with the usual calculation of gaussian fluctuations in the neighborhood of 
the superconducting transition leads to the expressions in Eq. _], i.e. Eq. [3] simply describes 
the dynamical fluctuations which lead to the Cooper instability as 7 — > 7» c . Here brackets 
( ) signifies the quantum expectation value. 

In the opposite limit, when the puddles are large with big values of the order parameter, 
one can neglect quantum fluctuations of the modulus of the order parameter and write the 
effective action in terms of fluctuations of the phase only 

s = £ G w/> j dTdT > \ _____ 1 + j drHj m] 

Hj[{<t>}\ = -(1/2) J2 Jij cos(0, - 0i) (4) 

The form of the effective action in Eq. |_is familiar from many earlier studies of the quantum 
dynamics of a system of superconducting grains linked by resistively shunted Josephson 



junctions (See, for example, Refs. _[ Q, __, 1121, [19|, |2_, l_ll-) In particular, the dynamical 



term proportional to in Eq. H]has the familiar Caldeira-Leggett |22j form and describes 

the quantum dynamics of the order parameter of an isolated puddle. In this case the origin 
of the dynamical term in Eq. HI is entirely different from that in Eq. [TJ it reflects the 
interaction of the phase fluctuations of the superconducting order parameter with quantum 
fluctuations of the electromagnetic field. In this case is the dimensionless effective 

conductance defined by injecting current into the i-th superconducting puddle embedded in 
the metallic host (ignoring the effect of other puddles), and measurin g th e voltage drop at 
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24j . associated with 



infinity. There is a further subtlety in 2D, as discussed in Refs. 
the fact that, G^ , as so defined, vanishes logarithmically with the size of the system. 
To handle this problem properly, one needs to consider corrections to the dynamical term 
in the effective action, Eq. HI When this is done, the result is equivalent to identifying 



G^ ~ V G<id- 

We will show below that in different situations either electromagnetic fluctuations or the 
Cooper instability can make the dominant contribution to the quantum dynamics of the 
order parameter . 



B. Josephson couplings between puddles 

We will see that near criticality the typical inter-puddle distance is larger than their size. 
The other generically important aspect of the problem is the dependence of the Josephson 
couplings between puddles on their separation, r-j — rj, which is long-ranged (power- law) in 
the limit that the temperature, T — > 0. Specifically, the coupling between small puddles in 
Eq. [1] is, up to logarithmic corrections which we will discuss later, of the form 



where v is the density of states in the normal metal, Lt = \/ D tr /T is the coherence length of 
normal metal which diverges as T — > 0, D tr is the electron diffusion coefficient, and Cjj(r^, r^) 
is, generally speaking, a complicated (random) dimensionless function of the coordinates. In 
the cases considered in Sees. II and IV, the value of C is determined by the average properties 
of the "normal" metal between puddles. In particular, in Sec. II, CV, is mostly positive, and 
so can be approximated by its average value, C^, while for the d-wave superconductor treated 
in Sec. IV, has a random sign, but at long distances, this sign is entirely determined by 
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J,lj = J(Yi,Yj) OC Cj 



'I] 



D 



cxp 



the character of the puddles, i and j, and is independent of the distance between puddles. 
By contrast, in the cases with a magnetic field considered in Sec. Ill , CV, is determined by 
the random quantum interference between different paths through the normal metal. As a 
consequence, has a random phase. 

In the limit of large puddles, the functional dependence of the Josephson coupling in 
Eq.H 



~ D tr 
J ij oc C%j 



R 



D 



R? |rj — ij 



D 



exp 



(6) 



on distance is the same as in Eq. [51 Here the random function has properties similar to 
Cij, and R is the typical size of the puddles. 



C. Susceptibility of an individual puddle 

The susceptibility of an individual puddle can be expressed in terms of the correlation 
function Xi(u;) (Eq. [3]) as Xi = Xi{ u — 0)- Its value depends on the puddle size. 

Let us start with the case when the puddle is large, so the dynamics of the order pa- 
rameter is determined by the effective action given by Eq. H]with Jij = 0. The implications 
of this effective action are best appreciated by interpreting imaginary time as a fictive spa- 
tial dimension, making the single puddle problem equivalent to the classical inverse-square 



XY model 



25 



26j at an effective temperature T e ^ = 1/G[ e ^\ The long-time correlation 
functions of the inverse X-Y model have been calculated in various ways, and are well un- 
derstood. The characteristic decay time depends exponentially on 1/T e ^ , and the dynamic 



correlation function has a power-law fall-off 271 ]. 



(A,'(0)4, W )HA/{ [^Pforr«r, (?) 

[Tj/r] for r > tj 

where X{ is a non-universal exponent Xi = T e ^ /(2n), the relaxation time, 

n ~ exp[ZG (e//) ] (8) 
and Z = 2tt 2 . The susceptibility is thus 

A exp[ZG 4 {e//) ]. (9) 
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However, in 2D, due to the dimension specific subtlety [23|, [24| discussed in Subsection IA, 
we must identify G^^' ~ VG2D, so that 



Xi 



A exp Z'^G^ (10) 



where Z' is another number of order 1. 

Let us now turn to the case when the modulus of the order parameter on a puddle is small 
and its dynamics are determined by Eq. [TJ There is a complicated, and for our purposes 
not terribly important crossover that occurs for puddles which are right on the verge of a 
mean-field transition, (7 — j ic \ -C 1. Ignoring such puddles, there are two distinct asymptotic 
behaviors that are readily deduced: 

. [aiilic-l)]' 1 for ( 7ic - 7) > 0, 

Xi ~ s „ I 11 ) 

Ai exp[ Z'% A* 2 ] for 7ic - 7 < 0, 



where, in this expression, |A$| = Aq-^/t — lie is the mean-field amplitude of the supercon- 
ducting order on puddle % and Z" is a renormalized relative of the same factor Z defined in 
Eq. [91 above. We now sketch the derivation of this result. 

The result for a puddle that does not support a mean-field solution, i.e. for -y ic > 7, is 
easily obtained by evaluating Eq. [3] at u = 0. In the opposite limit, (7 — 7ic ) > 0, there is a 
well developed mean field value of the order parameter on the puddle 

Aj = A OA /7-7 ic exp[z0 i ], (12) 

which has an arbitrary phase (pj. To the extent that modulous fluctuations can be ignored, 
this problem is precisely equivalent to the large puddle problem, with the role of G- 6 ^ 
played by PjA'j on each grain. Although large amplitude modulus fluctuations are relatively 
costly, because the resulting expression for the susceptibility depends exponentially on |Aj| 2 , 
they cannot be neglected. However, since the modulus of the order parameter appears 
exponentially in the expression for Xi, h is clear that the neglect of modulus fluctuations is 
not reasonable. However, they do not alter the asymptotic qualitatively, but rather result 
in a renormalization (reduction) of the factor Z" in Eq. [TH 

The most important feature of Eq. [11] is that the susceptibility increases exponentially 
as a function of 7 — 7j C and of the volume, Vi, of the puddle, % ~ exp[Z"z/A Vi(7 — 7i c )]- 
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D. Determination of the quantum critical point 



We now outline the procedure for determination of the location of the quantum critical 
point under these circumstances. 

Quantum fluctuations necessarily destroy the superconducting order in an isolated pud- 
dle. Thus, although the superconducting susceptibility of an individual puddle, \ii can > 
under some circumstances, be large, the transition to the globally phase coherent supercon- 
ducting state is ultimately triggered by the Josephson coupling between puddles. Let us 
introduce a dimensionless coupling between two puddles, i and j, 

Two puddles fluctuate essentially independently of each other if \X it j\ <C 1, and they are 
phase locked to each other if \Xij\ ^> 1. The transition to a globally phase coherent state 
occurs as a function of 7 at the critical value, 7 = j c , at which an infinite cluster of puddles 
is coupled together by links with JQj ~ 1. For an ordered array of p uddles, the quantum 



superconductor-metal transition was discussed in this light in 16l. l23l. |24J|. 



In disordered systems, the nature of the phase transitions described by the effective action 
in Eq. [1] depends on the distribution of the parameters, 7 ic , (3j, G- 6 ^ and J^, and these in 
turn are somewhat different in the various cases we treat below. 

However, what is common to the cases we will analyze is that, according to EqfjJ] the 
susceptibilities of the puddles depend exponentially on the parameters of the action Eq{TJ 
Thus in a generic situation in the neighborhood of the transition, the distribution of \i 
is extremely broad, and at criticality, rare puddles with exponentially large susceptibilities 
play a special role. In this case, the critical point can be identified by finding the set of 
"optimal puddles" which lie on the critical links of "the percolating cluster". This will be 
done in a way analogous to Mott's approach to the theory of the variable range hopping 
conductivity (See for example |28l|). 

Specifically, the optimal puddles are those in which j ic lies in an interval, ■y opt — ^ op t < 
lie < lo P t + ^lopt- Here both the optimal value, 7^, and the width of the interval, A7 opi , are 
determined by maximizing the quantity X opt = xlpt^opt with respect to these parameters, 
where Xopt is the susceptibility of a puddle with -y ic = 7 op4 , and J opt is the typical value of 
the Josephson coupling between two nearest-neighbor optimal puddles. Finally we find the 
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critical value of 7 = 7 C from the requirement that, after maximizing, max{X opt } ~ 1 . In 
the following sections we consider several examples of this program. 



II. A RANDOM MIXTURE OF METAL AND SUPERCONDUCTOR 



As a first case, we consider a random set of s-wave superconducting grains is embedded 
in a normal metal host with no magnetic field or magnetic impurities. We identify the 
superconducting grains as regions in which the effective interaction between two electrons 
in the Cooper channel is attractive (As > 0), while in the normal metal the interaction 
is repulsive (Ajv < 0). This system exhibits a metal-superconductor transition when the 
appropriate average value of A(r) changes sign, although the parameter kpl can still be 



arbitrarily large 



16 



24j . Some aspects of various 



closely relatec 
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24 



29 



prob 



30 



31 



ems have been 



321 ] . However, 



previously been analyzed using a variety of approaches j 161, 
we are able to obtain a much complete picture than has been obtained previously. Moreover, 
this problem serves as a useful warmup as it provides a simple explicit example of how 
the character of the optimal puddles and the nature of the quantum phase transition are 
determined from the present quantum percolation analysis. 

To be concrete, we will assume that the diameters of the grains, Ri, are random quantities 
characterized by the Gaussian distribution 



P{Ri 



N 



exp 



(Ri - Rf 
2a 2 R R 2 



(14) 



f 2ira R R 

where the average radius is R, the dimensionless variance <tr <C 1, and the total concentra- 
tion of grains is N. In the notation of the previous section, we can identify 74—7 = Ri/R c —1, 
where R c ~ £0 is the critical radius for the existence of a mean-field solution and £0 is zero 
temperature coherence length in the superconductor. 

Expressions for the susceptibilities of individual grains, in various limits, are given in Eqs. 
ElUDl andEJ 

The value of the Joshepson coupling between two superconducting grains embedded in a 
normal metal depends on whether the Andreev reflection on the N-S boundary is effective or 
not. When the puddles are larger than the coherence length determined by the magnitude 
of the order parameter in the puddle, one finds 



Jij ~ 



G 



D 



tr 



V 



eff 



R2 \r i -r j \ D [l + 2X N \\n 2 (\r i -r j \/R)\] 
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exp 



-T 



(15) 



by solving the Usadel equation. (See for example 33J). Here V is the volume of a grain, 
V ~ R D . In the opposite limit, when the value of the order parameter on the puddle is 
small and Andreev reflection is ineffective, the coupling can be computed from perturbation 
theory (See for example 



10| : 



vV 2 



Ti-Tj\ D [l+2|A JV |ln 2 (|r i -r J -|/ J R)] 



exp 



(16) 



Due to the s-wave character of the superconducting order and the fact that the system 
is time-reversal invariant, and so long as certain effects of strong correlations in the metal 
{jyj] can be ignored, is always real and positive. Moreover, if kpl ^> 1, the mesoscopic 
fluctuations of the magnitude of J„ are small compared to the average, and can thus be 
neglected. 

Our goal is to determine the critical concentration of grains N c at which the 
superconductor-metal transition occurs. We can identify various regimes depending on val- 
ues of the parameters <tr, (R — R c )/R c , G e ff, an d the Ginzburg parameter g = vR® A . (g is 
roughly the number of electrons within energy window Ao on an individual superconducting 
grain. Notice that, in the cases of interest here, where R c ~ £0, it is always the case that 
g ^> 1.) We analyze some representative cases: 

A) Monodispersed superconducting grains: Let us start with the case where <jr — > 
0, so all puddles are essentially the same size Ri = R. At T = 0, it is simple to see that: 

i) If (R c — R) < 0, individual puddles are not superconducting, even at mean-field level. In 
this case, which corresponds to the uniform "Cooper limit," the superconducting transition 
is, to first approximation, mean-field like and it occurs when the average interaction strength, 

X = NVX N - (1- NV)\X N \ (17) 

changes sign (from attractive to repulsive). Thus, the critical concentration is N C V ~ 
I'VwI/l/Vs + [ Aat I ] - This estimate neglects the spatial variations in the local concentration of 
superconducting grains; even when N is, on average, smaller than this mean-field critical 
value, there occur regions in which the concentration of grains exceeds this critical value. 
These regions act as the superconducting puddles of a new level of analysis, which gives 
results similar to those discussed below in the case of larger gr. 

ii) If 1 3> (R — R c )/R > 1/g, there is a well-developed mean-field order, A = 
^{R-R C )/RA on each grain, but the order parameter has a magnitude small compared 
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to Ao- In this case, Eq. [T] governs the dynamics, and we can make use of the expressions 
Eqs. [11] and [16] for x an d Jij, respectively. Since the grains are typically a distance of order 
N~ l l D apart, the dimensionless coupling between two neighboring grains is 

X ~ ^gR D N^(R - R c )/R c exp[Z"g(R - R c )/R c ] J . (18) 



This coupling is larger than 1 when N exceeds the critical density, 



gRP 



R — R c 

Rr 



1/2 



exp 



-Z"g 



R — R r 



Rr 



at T = 0. (19) 



iii) If (R c — R)/R > 1, the grains act, more or less, like pieces of bulk superconductor. Here, 
the dynamics of the quantum fluctuations are governed by electric field fluctuations as in 
Eq. HI and consequently the same analysis leads to 



, R- 2 exp(-Z'y/G^) in D = 2 at T = 
N c ~ { ^ V (20) 

R- 3 exp(-ZG eff ) in D = 3 at T = 0. 

Let us turn now to the temperature dependence of N C (T). Conventional arguments 
suggest that at low enough temperatures, and arbitrarily close to the zero temperature 



critical point, there is a universal quantum critical regime |35l. [36j, where, for example, 
(N C (T) — N c ) ~ N C T X , where x signifies an appropriate universal critical exponent. This 
regime, if it exists, applies only so long as 3> A^ c _1//Z) , and so is exponentially nar- 
row. Beyond the quantum critical regime, there is a broad range of temperatures in which 
N c l l D ^ Lt S> R, where the fluctuations are highly quantum mechanical in the sense that 
one can still neglect the T-dependencies of x(T). In this case there are two sources of the 
T-dependence of N C (T) : a) the classical fluctuations which destroy the coherence between 
puddles when Jy|Aj||Aj| ss T, and b) the T-dependencies of Jy(T), given by Eq. [T6] The 
relative importance of these two effects depends on the dimensionality of space D, and the 
value of the parameter g. The second mechanism dominates the T-dependence of N C (T) in 
the 2D case at arbitrary T, and in the 3D case in the wide interval of temperatures where 
Lt/R > g. Then the criterion X opt ~ 1 corresponds to a typical distance between puddles 
of order Lt, and hence 

N C (T) ~ — Lp, R « L T < N c (0y 1/D , (21) 

(In this article we will ignore relatively small temperature interval Lt > R where in 3D the 
temperature dependence N C (T) is determined by the first mechanism.) Notice, subtleties, 
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such as whether the superconducting state has long-range order (the 3D case) or only power- 
law order (the 2D case), do not affect the validity of this estimate. 

B) Optimal radius grains: If the variance a R is not too small, the susceptibilities of 
individual grains Xi have an exponentially broad distribution. As a result, at T = the 
transition point occurs when R < R c and is determined by relatively rare "optimal" puddles 
with anomalously large values of (R4 — R c )/R c , and consequently with exponentially large 
susceptibilities. However, as we shall see, if R is too much smaller than R c , the optimal 
grains become so rare that, yet again, a new regime occurs in which the optimal puddle is 
formed in regions with an anomalously large concentration of sub-critical grains. 

Let us focus on those grains with radii, Ri, within AR opt of a still to be determined 
optimal radius, R opt > R. (It can be shown that the relevant range is AR opt ~ u R R). We 
will ignore puddles which do not belong to this optimal set since puddles with much larger Rj 
are extraordinarily rare, and those with much smaller Rj have much smaller susceptibilities. 
Under the assumption (to which we will return below) that the optimal puddles are still 
small enough that Eq. [I] applies, the concentration of the optimal puddles is 

N opt ~ Na R exp[-(R opt - R) 2 /2a 2 R R 2 \ (22) 

so, according to Eqs. flTfTBl andl 



X opt ~ [N opt g fexp 
Maximizing Eq. [23] with respect to R opt gives 



// _ (Ropt Rc) 



Z"g 



Co 



(23) 



R^r ] = R+Z"gR 2 al/^ (24) 

and max{N opt } ~ A^exp [— (Z"gRcr R ) 2 /2^ <C A" (i.e. most puddles are smaller than R opt 
and hence play no direct role in the transition) . Finally, using the criteria X opt ~ 1 , we find 



Nc ^R^- R eXP 



Z"g ( Z"gR*o R 
tx c — It — 



at T = 0. (25) 



Co V 2£ 

We note that the critical concentration in Eq. [25] can be extremely small as a consequence 
of the fact that rare larger than average puddles contribute significantly to the global phase 
coherence. 

Let us now discuss the limits of applicability of Eq. [25j It is manifestly necessary that 
a R be small enough that N C R D 1, i.e. that 

(R c -R)> Z"gR 2 a 2 R /2^ . (26) 
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Note that that the same criterion leads to the inequality AR opt ~ Ran <C {Ropt — R), which 
justifies the saddle point approximation. 

The temperature dependence of N c can be obtained in similar ways as in A, above. 
When Lx 3> N^ D , the behavior is presumably governed by universal properties of the 
quantum critical point. However, for N^ D > L T ^> R, the temperature dependence of 
N c again derives from the temperature dependence of Jjj, i.e. that the Josephson coupling 
falls rapidly to zero at distances large compared to Ly. Here, the character of the optimal 
puddles is still determined by Eqs. [23] and EH but the critical concentration is determined 
by the condition that the distance between optimal puddles is typically of order Lt, i.e. 



The resulting phase diagram is shown schematically in Fig. 1. 

C) Large puddles: The expressions in B were derived under the assumption that the 
optimal grains are sufficiently small that the mean-field order parameter has a magnitude 
small compared to Ao, and consequently that the susceptibility grows exponentially with 
the radius of the grain, according to Eqs. [2] and [TTJ However, if or is sufficiently large, 
the optimal grains get to be large enough that R opt — R c > A£ (where A is of order 1), 
and consequently |Aj| ~ A . In this limit, the quantum dynamics of the order parameter is 
determined by quantum fluctuations of the electromagnetic field, and Xi is determined by 
the normal state conductance of the metal, G e ff, as in Eqs. l9l and [TU1 We can estimate 
the conditions for this by extrapolating Eq. [21] to the point R op t — R c = £o ; from which we 
deduce that the optimal puddles are "large" when <j\ > (£ / R) 2 : (1 / Z" 'g) . 

In this limit, since \ grows with size of the grain only relatively weakly (log[x] ~ R D ~ 2 ), 
while the density of grains of large size falls exponentially with their volume, the density of 
optimal puddles is simply the density of grains with radius larger than R = R c + A£ , 



With N[ arge playing the role of N opt , and with the expressions in Eq. [9] and [10] for the 
susceptibility, the same analysis can be applied as in part B to obtain N c . For instance, in 



N C (T) ~ LTjP exp (Z"n c Ra R ) 2 /2H 



(27) 




(28) 



3D 



N c ~ R~ 3 exp[-ZG eff + (R— R c f /2R 2 a 2 R ]. 



(29) 
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FIG. 1: Schematic phase diagram for the case considered in Section IIB, in which a concentration 
N of superconducting grains with a distribution of sizes is embedded in a metallic host. The 
solid line represents a phase transition. The dashed line represents a crossover where a "small 
fraction" of the sample first supports a mean-field solution - although near this line, global phase 
coherence is destroyed by quantum and thermal phase fluctuations, measurable manifestations 
of local superconductivity onset below this line. The dotted line represents the expected phase 
boundary in the Cooper limit, where mesoscopic fluctuations are ignored and a single uniform 
effective interaction between electrons is assumed as in Eq. [T71 

We would like to stress that Eq. [23 holds only at exponentially small temperatures for which 
Lt ^> N c l l D . In the opposite limit, puddles with Ri > R D are irrelevant, and the Cooper 
instability contribution represented by EqJTT] dominates the physics of the phase transition, 
even for relatively large values of a R . 

III. S-WAVE SUPERCONDUCTOR IN A MAGNETIC FIELD 

A. A magnetic field perpendicular to a superconducting film 

A magnetic field, H, applied perpendicular to a metallic film (D = 2) couples primarily 
to the electron's orbital motion. In this case the Zeeman coupling of the magnetic field to 
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the electron spin can be neglected. To the extent that mesoscopic fluctuations of the order 
parameter can be neglected, the problem of s-wave superconductivity in a magnetic field 
was solved by Abrikosov and Gorkov. In this approximation, at H^ = $ /7r^, the order 
parameter can be represented as a superposition of wave-functions in the first Landau level, 

<f>(r) = -j=t= exv(-r 2 /L 2 H ), (30) 



where $ G = hc/2e is the flux quantum, L H = y ^ /2ttH^ is the magnetic length. (We 
consider the "dirty limit" £ 3> I 3> kp 1 where £ = y^D tr / A ). Roughly speaking, the 
same form of the wave-function applies even when mesoscopic fluctuations are taken into 
account. This simplifies the analysis in that it implies that, near the point of the transition, 
the puddles have a typical size, Lj ps Lh- However, the critical magnetic field Hi varies 
randomly as a function of position, so that, in the notation of Section I, we can identify 

_ (Hi ~ H) 

7<-7 = — -<5j-. (31) 
We will assume that the distribution of Hi is approximately gaussian 



1 , (Hi-H c 



|2 



^7^ exp[ -^pf' (32) 



and is characterizing by the average H c and a dimensionless variance cr#. (This ignores the 
existence of long, but for present purposes irrelevant tails of the distribution produced by 



mesoscopic effects 37].) We assume that on <C 1, and thus that H^ ps H c . Generally, 
there are two contributions to the variance an'- one contribution comes from classical fluc- 
tuations in the strength of the local scattering potential and one from non-local quantum 
"interference" effects, 

a H = a {int) + a (d) . (33) 
The classical contribution is due to random fluctuations of the concentration of impurities, 

a^(L) ~ {A/ L S ) (34) 

where A <C Lg is the correlation length of the disorder potential. The electron interference 



contribution is 



38j: 



o-M) _ \jG %D < 1. (35) 
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Note that, although the interference term is independent of puddle size, and hence is the 
larger term for big enough puddles, for large G^d there is a parametrically wide range of 
puddle sizes for which the simple statistical variations in impurity concentrations dominates 
the variance of local critical fields. 

The configuration dependent (mesoscopic) variations in the Josephson coupling, J^ are 
more important here than in the previous example. One can see this by noticing that at 
large |r$ — rj| (up to possible logarithmic corrections) the averages 



Jij ~ 



uL 2 H d- 



T 2 



exp 



Jij ~ 



G 



D 



tr 



eff] 



exp 



|r,- — r 



|r,- — r 



(36) 



are much smaller than the variances 

1/2 



I 7--I 2 



I 7--I 2 



1/2 



J 2 

Vpl 
D tr 



Hi 



exp 



r, ; - r 5 



exp 



(37) 



where O indicates the average of O over realizations of the random scattering potential., 
and vf is the Fermi velocity. As a consequence, at long distances the Josephson coupling in 
Eqs. [T]and|l]is of the form 38] 



Jij ~ 



Jij ~ 



FijuL 2 H d- 



T 2 
H 



r, - r, 



exp 



L 



T 



F 



D 



tr 



exp 



L 



T 



(38) 



where = Fij^rj), and F i3 



domly in phase, and |Fy 



Fij(ri,Tj) are dimensionless functions which varies ran- 
~ 1. 



To find the critical magnetic field H c , we employ the same optimization procedure we 
used in the Section 2. We introduce an interval in the space of Hi which is centered at H opt 
with width AH ~ <JhH c . We will see that for sufficiently large values of an, the distance 
between the "optimal puddles" is large enough that the Joshepson coupling between puddles 
is dominated by the mesoscopic contribution given by Eq. [3S1 

Depending on the value of G2D = e 2 b , D tr d, the quantum dynamics of the order pa- 
rameter is determined by either the Cooper instability or the quantum fluctuations of the 
electromagnetic field. 
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If (Hopt — H)/!!^ -C 1, and hence A opt <C A , the quantum dynamics of the order 
parameter is determined by the Cooper instability, so we can use Eqs. 2 and 5 to obtain 



X opt ~ a H exp 



rrll 

Z 9 



H pt — H \ {Hopt — H c ) 



\2 



H r . J o\El 



at T = 0. (39) 



Here a# ~ A /E F , and g = vA^d ~ G 2 d- As before, H op t is the value which maximizes 
this expression, and the true critical field is then determined as the value of H = H c at 
which X op t = 1: 

Hopt = H C [1 + {ZG 2D /2)a 2 H } ; H c = H c [l + (ZG 2D /4)4] . (40) 

These expressions are self-consistent so long as 

2 2 

\fZG~2D ZG2D 

where the first inequality guarantees that {H c — H c )/H c 1 and the second that th e op timal 



puddles are dilute. The expression for H c in Eq. HQ] has been obtained in Ref. |l8j] by a 
different method. 

In any puddle for which {Hi — H)/H^ ^ 1, the superconducting order is "fully devel- 
oped," Aj w Ao, so the dynamics of the order parameter is determined by the quantum 
fluctuations of the electromagnetic field. In this case, the susceptibility of the puddle de- 
pends only on G 2 d as in Eq. [TO], and is independent of Hi. Since the probability of finding 
such a puddle decreases with increasing Hi, the optimal puddles of this sort are those with 
Hi w H op t ~ H. The corresponding dimensionless coupling between these puddles is thus 

{H op t — H<y 2 



X pt ~ exp 



2Z'^/G 



2D 



(42) 



These puddles are always dilute compared to their size so long as [H — H c ] > ohH c . The 
critical value of H, determined by the condition X opt ~ 1, is 

H c = H C {1 + V2a H Gl^) at T = 0. (43) 

The issue of whether the global superconducting properties are dominated by weakly 
superconducting or "fully developed" puddles is settled by determining which of the expres- 
sions for X Q pt in Eqs. |39] and H2] give the largest value. In particular, the critical field is 
determined by the larger of the values given by Eqs. HQ] and H3] 
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Since in both these cases, the puddles are dilute, the phase of the Josephson couplings 
between optimal puddles is random. Thus, the T = ordered state is glassy, and is an 
example of a "gauge-glass." While this phase does not have long-range order, it supports a 



non-zero Edwards- Anderson order and is generally thought to have zero resistance. [411. 142] 
In 2D, however, there is no ordered state at non-zero T, so there are only crossovers as 
a function of H and T. There is a characteristic field H*(T), such that for H < H*, the 
coupling between optimal puddles is large (in magnitude) compared to the temperature; 
here, the resistivity is, presumably, due to some form of variable-range-hopping of vortices, 
and so decreases exponentially with decreasing T. Clearly, H*(T) — > H c as T — > 0. It is only 
weakly T dependent at low T, but in the temperature range such that is smaller than 
the typical spacing between optimal puddles, but large compared to the puddle diameter, 
H* is determined by the condition that the concentration of superconducting puddles must 
exceed L^ 2 ; this leads to the unusual T dependence: 

H*(T) ~H e \l + 2a H ln 1/2 (L T /L H )] . (44) 

The resulting schematic phase diagram is shown in Fig.2. Note that there are a series of 
additional crossovers that we have not addressed here, and which are not shown in the figure, 
which occur at fields greater than H*. These crossovers characterize various energy scales in 
the anomalous metallic phase proximate to the superconducting glass. More of the physics 
of the anomalous metal will be addressed in future studies. 



B. The case of a parallel magnetic field 

We now consider the opposite limit, in which the coupling of an applied magnetic field 
to the electron spin (Zeeman coupling) is significant, and the coupling to the orbital motion 
of the electrons can be neglected. To a good approximation, this can be realized in a thin 
film of an s-superconducting metal in which the superconductivity is destroyed by an in- 
plane magnetic field H\\ . The situation in this case critically depends on the value of the 
parameter Aqt so (or, in other words, on the atomic weight of the metal) where l/r so is the 
spin-orbit scattering rate. In the case of relatively strong spin-orbit coupling, A t so <C 1, 
on mean field level and in the absence of mesocopic fluctuations, the transition is second 
order. In this case, the effect of mesoscopic fluctuations on the character of the transition is 
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qualitatively similar to the transition in the perpendicular magnetic field considered in the 
previous subsection. In the case A r so <C 1, however, the situation is very different because 
on the mean field level the transition is first order. In this section, we will consider this case, 
and to simplify the discussion, we will consider it in the limit of zero spin-orbit coupling, 
AoT so — > oo. In this limit, the microscopic physics responsible for the emergence of a puddle 
state is quite different, and in particular the energy associated with the formation of puddles 
is larger than in the previous examples. Therefore the various manifestations of the physics, 
and especially the glassy character of the phase at intermediate fields, is more robust. 

If mesoscopic fluctuations are ignored, the zero temperature transition is first order, with 
a discontinuous jump in the spin-magnetization density from m = (in the superconducting 
state) to m = XspH\\ in the metallic state, where Xsp — vh^i 2 B is the normal-state spin (Pauli) 
susceptibility. The critical field is given by the well-known Chandrasekar-Clogston limit, 

< = A /hfi B . (45) 

(In the disorder free case, there might appear a narrow regime of fields near H^) in which 
a partially polarized Fulde-Ferrel-Larkin-Ovchinikov (FFLO) state occurs, but this is not 
relevant in the case / <C £o, considered here.) 

In disordered systems the critical magnetic field H c n(r) exhibits spatial fluctuations. In 
the absence of spin-orbit scattering the value of H$ given by Eq. 05] is independent of I, 
which means that there is no classical contribution to the dimensionless variance cry, which 
is, instead, determined entirely by mesoscopic interference effects. Thus, the dimensionless 



variance is 



39] 



(H c \\ — H c \\) 



1/2 



H c w G 



' (46) 



2D 



and the correlation length of H c \\(r) is of order £ - 

According to general theorems j43|, in 2D quenched disorder destroys first order tran- 
sitions. Rather, near a putative first order transition, a domain structure occurs with a 
characteristic domain size Ldom which is exponentially large in the small disorder limit, 
Ldom ~ ex P[Z"0-/°'\\) Z ]- I n some cases, the putative first order transition is simply smeared 
and replaced by a crossover which becomes increasingly sharp as the disorder becomes 
weaker. However, in the present case, since there is clearly a superconducting phase for 
small enough H\\ and (as we shall confirm) a non-superconducting phase for large enough 
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Hh, there must still be a sharp, continuous quantum phase transition at a shifted critical 
field, H 4 . 

Specifically, in the present case, the superconducting domains are regions with magneti- 
zation m near 0, and with the local magnitude of the order-parameter, |Aj| ~ A , while the 
metallic regions have m ~ XspH\\ and miniscule magnitude of the superconducting order. 
The volume fraction of the two phases is a function of H\\; it is roughly a 50-50 mixture 
when H\\ ~ H^), and the superconducting fraction decreases monotonically with increasing 
H\\. However, global phase coherence is not lost at H\\ ~ H^< where on the mean field 
level the superconducting fraction first fails to percolate. Rather, as in the other problems 
we have examined, it occurs when the Josephson coupling between superconducting regions 
becomes sufficiently weak, which in turn occurs when the superconducting fraction is small 
and the superconducting regions far separated. 

Because the superconducting regions have a characteristic size large compared to £ , 
and the magnitude of the order parameter is large, the dynamics of phase fluctuations 
is determined by electric field fluctuations, and consequently (according to Eq. [10]) Xi ~ 
Aoexp [Z'y/G^\. 

To determine the distribution of Josephson couplings, we note that in an SNS junction, 
when the normal part of the junction is partially spin polarized, 40j the Josephson coupling 
oscillates in sign as a function the coordinates. Specifically, at T = 0, 



■r. ,, G 2 DD tr , r -i, 

J(r, r) ~ — — exp( ) cos 
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J(r, r' 



r-r'|" L H \\ L H \\ 
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tr 



elr — r'l 2 



(47) 



J(r, r') ~ F(r, r')- -J cos( — 

r — r r L 



where Ljj\\ = \jD tT / /ZHy, and F(r,r') is a sample specific function (\F\ ~ 1) which has 
random variations both in modulus, and in sign. 

The mesoscopic fluctuations of J again dominate the average at distances large compared 
to Ljj». Thus we can estimate the critical magnetic field H c u at which the zero temperature 
phase transition to the metallic phase takes place. For H» > H c n , the probability of finding a 
superconducting puddle is ~ exp — [H^ — H c \^ 2 /2a?<H\ . As a result, following the same 
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FIG. 2: Schematic phase diagram for the cases considered in Sections IIIA and IIIB in which super- 
conductivity in a thin film of an s-wave superconductor is destroyed, respectively, by application of 
a perpendicular or a parallel magnetic field. The crossover scale, H*, (indicated by the dot-dashed 
line) is described in the text. The "gauge glass" refers to a zero temperature phase in which the 
resistance vanishes as T — > 0, but for which there is no finite temperature phase transition. The 
mean-field phase boundary has a continuous portion (indicated by the narrower dotted line) and a 
first order portion (indicated by the heavier line) separated by a tricritical point (indicated by the 
solid circle). 

line of reasoning as in previous sections, 



Notice that the average Josephson coupling in Eq. 271 itself, oscillates in sign, so even 
when the superconducting puddles are closely spaced, the Josephson couplings are random 
in sign. The resulting glassiness of the superconducting state is more robust than in the case 
of the perpendicular magnetic field because of the large size of the superconducting domains, 
the fact that the magnitude of the order parameter is fully developed within each domain, 
and due to the fact that the randomness in sign is not solely a long distance subtlety. In 
the absence of spin-orbit coupling, the magnetic field does not couple directly to any orbital 




H r \\ — H, 



(48) 
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degrees of freedom, and hence the glass phase can be precisely characterized as in an XY 
spin-glass, in which the ordered state supports spontaneously generated equilibrium orbital 
currents. 

In the presence of spin-orbit coupling, a sharp definition of the glass phase based on 
the presence of orbital currents is not possible; even a normal disordered metal phase will 
support small orbital currents under these circumstances. As in the previous section, the 
superconducting coherent state near the critical field is some form of a gauge glass. 



IV. DESTRUCTION OF D-WAVE SUPERCONDUCTIVITY BY DISORDER 



We shall now consider the case in which, in the zero disorder limit, the superconducting 
state is a BCS state with d-wave symmetry due to a weak attractive interaction in the 
d-wave particle-particle channel. 

In the d-wave case, it is necessary to explicitly treat the dependence of the supercon- 
ducting order parameter on the relative coordinate. Specifically, in the absence of disorder 
and in a bulk sample, A(r, r') = A^/cjj (r — r'), where k^\t) changes sign under rotation 
by 7r/2, and is a short-ranged function, with range, b <C determined by the range of the 
effective attractive interaction. 

Although the "d-wave" notation is inherited from spectroscopic notation for an "Z = 2" 
irreducible representation of the rotation group, in a crystal, it refers to an appropriate 
irreducible representation of the point group. We will treat the case in which the point 
group has at least two distinct even parity one dimensional representations - a trivial one 
and a non-trivial one. For instance, in a tetragonal crystal, in addition to the trivial (s-wave) 
representation, there are three other even parity irreducible representations: a d x 2_ y 2-wave, 
a d XJ/ -wave, and a g-wave (which transforms like (x 2 —y 2 )xy). We consider the case in which 
in the zero disorder limit, there is an effective attraction only for one of these representations, 
which we will call simply the "d-wave." 

The most clear-cut manifestation of the d-wave nature of the ground state order parameter 



comes from "phase sensitive" measurements 
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45] of the symmetry of the order parameter. 



Specifically, in a corner SQUID 45] of the sort described in Fig. 4, in which the external 
circle is a conventional s-wave superconducting wire, the ground state of the system will 
contain a half-flux quantum trapped in the SQUID for the case in which the sample is a 
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d-wave superconductor, and no flux if it is an s-wave superconductor. 

The fact that, in the absence of disorder, A(r, r') changes sign under rotation makes 
the system much more sensitive to the strength of the disorder than a conventional s-wave 
superconductor. We will characterize the disorder strength by the electron mean free path 
I in the normal metal. What happens to the system in the presence of relatively strong 
disorder depends on the sign of the interaction in s-channel. If the interaction in the s-wave 
channel is attractive but much weaker than the attraction in the d-wave channel, then when 
the disorder is weak enough (/ > £ ), the d-wave state dominates, but when the disorder 
strength increases enough to destroy the d-wave superconductivity (I < £o), the system 
undergoes a phase transition to an s-wave state (See for example 46j). The s-wave state is 
destroyed only when kpl ~ 1. However, in this article we consider the more interesting case 
in which the interaction in the s-channel is repulsive, so when disorder suppresses d-wave 
superconductivity, it drives the system to a normal state when the mean free path is still 
relatively large, kpl 3> 1. This case may be relevant, for instance, to the destruction of 
superconductivity in the "overdoped" high temperature superconductors. 

In the (conventional) approximation in which spacial fluctuations of the electron mean 
free path are neglected, d-wave superconductivity is destroyed when I ~ l — 1.78£ Q . Thus 
disordered d-wave superconductors are another example of a system which may have a 
quantum superconductor-metal transition in a situation in which the conductance is large. 
This case exhibits both similarities and differences with the cases we have already considered 
in Sections 2 and 3. 

In the presence of disorder, a material has no particular spatial symmetry, and so the order 
parameter cannot be said precisely to have any particular symmetry at all. Nevertheless, in 
bulk samples, symmetry is restored upon configuration averaging. It is therefore legitimate 
to ask questions concerning the global symmetry of the order parameter. Hence, we can ask 



whether A(r,r'), or J-(r, r') have d-wave or s-wave symmetry. Here the overline stands for 
a configurational averaging, and !F(r, r') = !F(r, r', t = t') is the anomalous Green function 
which is connected to A(r, r') by the interaction constant. 

It is important to realize that it is possible (indeed, as we shall see, inevitable near 
criticality) to have a situation in which the local pairing is "d-wave like" and yet the global 
superconductivity has s-wave symmetry. In fact we will show that there are at least two 
quantum phase transition as disorder increases: the transition from d-wave to s-wave global 
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symmetry, and subsequent transition from s-wave superconductor to the normal metal. The 
existence of the second (d-s) transition is the main difference with the cases considered in 
previous sections. (In fact, we consider it likely that rather than a sharp d to s transition, 
there is an intermediate glass phase in which time reversal symmetry is broken and s and 
d-wave ordering coexist. However, we have not fully explored this scenario.) 

We propose several different definitions of the global symmetry of the order parameter: 
a) The best operational definition is provided by the result of a phase sensitive experiment, 
such as the corner SQUID experiment shown in Fig. 4. b) The quantity A(r, r') can be 
characterized as having d-wave or s-wave symmetry. It can also provide a definition of a 
state with coexisting order if it has mixed symmetry, c) A specific diagnostic for a globally 
s-wave component of the order parameter can be defined in terms of the local component of 
the anomalous Green function jF(r = r') = J-^(r). If we define P± to be the volume fraction 
of a sample where F^- S \r) has a positive or negative sign, respectively, then the system has 
an s-wave component if (P + — P~) ^ 0. These definitions are not equivalent under all 
circumstances. However, for the purposes of this article, we are not primarily interested in 
the most general definition of the global symmetry of the superconducting state. For the 
most part, we will deal with the interval of parameters in which all these definitions are 
approximately interchangeable. 

1. The d-wave to s-wave transition as a function of disorder 

This transition takes place in the region of concentrations where quantum fluctuations 
of the order parameter can be neglected, and therefore it can be understood on the mean 
filed level. As a warmup exercise, consider a cartoon picture of a system of superconducting 
puddles of a size large compared to £ an d of a rectangular shape which are embedded in a 
noninteracting diffusive normal metal (See Fig. 5). The rectangles are identical, and they 
are oriented either vertically, or horizontally in a random fashion. The order parameter 
inside the rectangles has d-wave symmetry, and the orientation of the gap nodes is assumed 
to be pinned by the crystalline anisotropy. 

In a d-wave superconductor, in addition to an overall phase of the order parameter, 
there is an arbitrary sign associated with the internal structure of the pair wavefunction. 
Specifically, we adopt a uniform phase convention such that when the phase of the order 
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FIG. 3: Schematic phase diagram for the case considered in Section IV in which a BCS (weak- 
coupling) d-wave superconducting state (in 3D) is destroyed as a function of increasing disorder 
strength. The dotted line represents a transition between d-wave superconuctor and normal metal 
calculated by the conventional theory. In the presence of disorder, the labels "s-wave" and "d-wave" 
refer to the behavior of the system in a macroscopic phase sensitive measurement, as described 
in the text. (The negative slope of the boundary which divides globally d- and s-wave global 
superconductors shown in the figure can be justified only in the case when the electron interaction 
in the s-channel is attractive, so that the entropy of the s-wave superconductor is smaller than the 
entropy of d-wave. More generally, the slope of this boundary is to be determined.) 

parameter <pi = 0, this implies A(r, r') in puddle % is real and has its positive lobes along 
the (appropriately denned) y axis and its negative lobes along the x axis. 

It is obvious that at a high concentration of puddles, the order parameter in the ground 
state has global d-wave symmetry (See Fig. 5a.). However at small puddle concentrations, 
the situation is different. If the distances between puddles | — | R are much larger than 
the characteristic size of the puddles, R, the Josephson coupling between puddles inevitably 
favors globally s-wave superconductivity, even though the order parameter on each puddle 
looks locally d-wave like. In this case the mean field exchange energy of the system has a 
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form 



(49) 



where rji = ±1 are random numbers such that rji = 1 for a rectangle oriented in the x- 
direction and rji = — 1 for a y-directed rectangle. 

Eq. H9] represents the Mattis model which is well known in the theory of spin glasses 
[47I • The ground state of this model corresponds to 



Thus the distribution of exp(0j) between puddles looks completely random, as shown in Fig. 
5b. However the system is not a glass because it's ground state has a hidden symmetry, 
which in the present problem corresponds to to a global s-symmetry of the order parameter 
according to any of our proposed definitions! 

A qualitative explanation of this fact is as follows: The inter-puddle Josephson coupling 
originates from the proximity effect in the normal metal, which is characterized by the 
anomalous Green function !F(r, r'). Due to lack of symmetry at the boundary of a puddle, 
an s-wave component jF(r = r') = !F^ S > (r) 7^ of the anomalous Green function is generated 
in the neighboring metal. Specifically, at a normal metal-superconductor boundary, the 
sign of jF( s )(r) is determined by the sign of the d-wave order parameter in the k-direction 
perpendicular to the boundary (See Fig. 5). (Thus the sign of jF^^r) changes along the 
boundary of a puddle.) On distances from the boundary larger than I, the anomalous 
Green function becomes isotropic. In other words, only the s-component F^ s '(r) survives 
elastic scattering. It is this component which penetrates through the metal and carries the 
Josephson current between puddles. At distances larger than the size of the puddle (but 
smaller than |r^ — |) the quantity jF s (r) has a definite sign which is determined by an 
integral around the surface, which sign gives us the value of rji. 

On intermediate distances, the situation is more complicated. Areas with different signs 
of J-^ s '{t) mix in a random fashion. We argue that the most important aspects of this 
complex situation can be modelled by the following effective Josephson energy: 



i+3 

where J^- characterizes the strength of the exchange interaction between the d-wave compo- 
nents of the order parameter. Typically at small |r$ — jff > but at large |r, ; — the 



cos(^) = -Tji. 



(50) 




(51) 
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?IG. 4: A schematic picture of a phase sensitive "corner SQUID" experiment, introduced in Ref. 



451 ] . If the square piece of superconductor has global d-wave superconductor symmetry, then there 
is a magnetic flux trapped in the ground state of the system. Pluses and minuses inside rosettes 
indicate the sign of A(k) as a function of the direction of k 

coupling strength J^' decays more slowly than jff . Thus between the asymptotic d-wave 
dominated regime where the puddles are dense, and the s-wave dominated regime where 
they are dilute, it is likely that there is an intermediate region in which the system will 
exhibit spin glass features and/or coexistence of d-wave and s-wave ordering. In this article, 
however, we will not further explore this fascinating but complex aspect of this problem. 

While the above discussion was based on a cartoon model with regularly shaped puddles, 
we would like to stress that our conclusions do not rely on this. In particular, as it is 
qualitatively illustrated in Fig. 5c, Eq. H9] holds at arbitrary shape of the puddles provided 
that the typical distance between them is larger than their characteristic size. 

To quantify our conclusions, we compute the Josephson coupling between a pair of far 
separated puddles in two extreme limits, large puddles and small puddles: 

If the size of the puddles is large enough, the Josephson coupling has to be obtained from 



the solution of the Usadel equations (See for example [33() for the configuration averaged 



29 




FIG. 5: A qualitative illustration of the global d-wave to s-wave transition. Solid lines represent 
boundaries of d-wave superconducting puddles embedded into a normal metal. Pluses and minuses 
indicate the areas were the s-wave component of the anomalous Green function F s (r, r) is positive 
and negative respectively, a) The "cartoon" case in which the concentration of regular rectangular 
d-wave puddles is large, so the system has d-wave global symmetry, b) The case when the con- 
centration of d-wave puddles is small so the system has s-wave global symmetry, c) The case in 
which the concentration of d-wave puddles is small and they have arbitrary shapes. Here, they are 
shown embedded in a normal metal, and the system has global s-wave symmetry. 



anomalous Green function (r) = — 2sin#(e,r) in the metal, 

-^0(e,r) +iesin0(e,r) = 0. (52) 

Here jF^^r) is the Fourier transform of J-( s \r, t — t'), and the symbol indicates averaging 
over the random scattering potential between the puddles at given shape of the puddles. The 
boundary conditions for Eq. [52] at the normal-superconductor surface have been derived in 
Ref. [481 1 . They are valid as long as the size of the puddles is large and the Andreev reflection 
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on the puddles is effective. 

Since the relevant energy for computing the Josephson coupling, e ~ D tr /\v — r ^ | 2 , is 
much smaller than the value of the order parameter in the puddles, the boundary condition 
for 9(r, e) is independent of e, and depends only on the angle between the unit vector parallel 
to the direction of a gap node, and a unit vector, n(r), normal to the boundary at point 
r at the surface, : 

B (e, r) = f[a(r)], sin[a(r)] = n(r) ■ n A . (53) 

Here f{o) is a smooth, approximately odd and periodic function, f(a) ~ —/(—a), f{ot) ~ 
f(a + 7i"), which grows from f(a) ~ at a = 0, to /(a) ~ ±C f° r a = 7r /4- Here £ ~ 1. 

Solving Eq. [52] with the boundary conditions Eq. [53] and using the standard procedure 
of calculation of the Josephson energy we get 

i]i = sign | Jdsf (a) J (54) 

~(s) 

where the integral is taken over the surface of the ith puddle. Moreover, the value of 
in Eq. [49] turns out to be of the same order as in Eq. [4] It is this long range nature of the 
decay which ensures the existence of the phase in which the puddles separates by a large 
distance and the system has global s-wave phase coherence. 

In the second case, the Andreev reflection is ineffective and the interactions between 
puddles can be computed in perturbation theory. Thus we can write an analog of Eq. 



^ 8 --E^[ J ? A * A i+ c ' c -] ( 55 ) 

i¥=3 



where, Aj is (up to a sign, 77'- = ±1) the average of the order parameter over puddle 

, drdr' 



1, 

Ai = Vi I A(r,r')^, (56) 

J puddle i "i 

and rj[ is a random variable that we have introduced in Eq. [55] (and then cancelled in the 
definition of Aj). The strength of the Josephson interaction between s-component of the 
order parameter in puddles characterized by is (up to a numerical factor smaller than 
one) of order of Eq. [TBI Again, we have neglected in Eq. E51 the interactions between the 
d-wave components of the order parameter, since they fall faster with separation between 
puddles. Notice that in the fine tuned case of a fully symmetric puddle, Aj = due to the 
d-wave symmetry. 
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The important point is that, in both the small and large puddle limits, Eqs. [55j and H9l 
yield the same qualitative picture: at large inter-puddle distances the Josephson coupling 
favors s-wave symmetry. It now remains to show that, near the point of quantum SMT, the 
distance between optimal puddles is indeed much larger than their size. 

2. Globally s-wave superconductor to metal transition 

The quantum transition between a globally s-wave superconducting state and the metal 
does not differ in a crucial way from the transition which has been considered in Section 2. 
For reasons that should by now be familiar, near the critical value of the disorder the spacial 
dependence of the order parameter can be visualized as defining a system of far separated 
superconducting puddles with anomalously large values of the order parameter separated 
by large areas of the normal metal. In particular, this results in a smaller value of l c (larger 
critical magnitude of the disorder strength) for the destruction of superconductivity than l co 
which is given by the conventional theory. The difference between the present problem and 
that treated in Section 2 is that, to identify a set of optimal puddles, one has to characterize 
them by two generally independent parameters: the size of the puddles and the value of the 
s-component of the order parameter associated with an optimal puddle. 

As we have seen in previous sections, depending on whether the Andreev reflection from 
the superconductor-metal boundary is effective or not, two scenarios are possible. In the first 
case the amplitude of the Josephson energy is independent of the value of the order parameter 
in the puddles, and is given by Eq. [15j In the second case the amplitude of Josephson energy 
is given by Eq. [161 The susceptibility of an isolated puddle does not depend in any essential 
way on the symmetry of the superconducting order parameter, and as such, the analysis 
follows along identical lines as for the case of s-wave puddles, discussed in Sec. II. Thus, 
for small puddles, the susceptibility is determined by the Cooper instability Eq. [HI while 
puddles of radius R which is large compared to the coherence length, the susceptibility is 
determined, as in Eq. [9j by the effective conductance, G e ^ ~ R D ~ 2 , where R is the radius 
of the puddle. 

Let us consider a situation where the mean free path Z(r) is a random quantity which 
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exhibits classical spatial fluctuations with a Gaussian distribution 



= vkf xp[J W ] (57) 



characterized by the average I, a dimensionless variance 07, and a correlation length, A. 

To be concrete we consider the 3D case, and neglect the mesoscopic fluctuations of the 
mean free path of an interference nature. We also assume that the conductance of the metal 
is isotropic. 

The size of the optimal puddles is readily seen to be of order the coherence length R opt ~ 
iopt ~ i ~ £o\Av 'Q-opt ~ lo) j an d therefore the susceptibility of the puddles is given by 
the large puddle result in Eq. [91 (The self-consistency of this assumption can be checked 
starting with the assumption that the small-puddle expression in Eq. [TTJ can be used, 
and then determining the optimal puddle size - this procedure leads to the inconsistent 
conclusion that the optimal puddles are arbitrarily large.) 

The effective conductance that determines the susceptibility is the conductance of a region 
of characteristic linear dimension ^ opt \ this, in turn, depends on the local value of the mean- 
free path as G opt = G^ \/Iq/ (l opt — loj, where G$ is a conductance of a region of size £ - The 
probability to find a puddle of this size is determined by the variance of the mean free path 
averaged over it's volume, which is of order <7;(A/£) 3/ ' 2 . Thus we have 



X opt ~ exp 



(58) 



As usual, l opt is the value which maximizes Eq. [58l and the critical disorder, I = l c , is 
obtained by equating the result to unity. Although l c > Iq, as long as 07 is sufficiently small, 
(^0 — Icl/l^. 1- In this limit, the result of this procedure can be expressed as 

(l -l c )/lo~2ZGt a?(^) (59) 

Near the critical point I = l c , the distance between the optimal puddles 
Ropt exp[G^ 2 (£ /A) 3 / 2 /a7] ^> R opt is exponentially large. The assumption that the distance 
between optimal puddles is large which we made in the previous subsection is, thus, justified 
for / near l c . 

The temperature dependence of the critical value of the mean-free path, l c (T), can be 
obtained from similar considerations to those used to determine N C (T) in Section 2. 
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The resulting phase diagram for a d-wave superconductor in the presence of quenched 
disorder is shown schematically in Fig. 3. At large disorder (I < l c ) the system is in the 
normal metal phase. At Z > L the system in a state with a dominant s-wave component 
of the order parameter, and a d-wave component, whose sign is locally slaved to the s-wave 
component, in a way which varies randomly in space. At still larger values of I > Iq , there 
is a dominantly d-wave state, in which the s-component has a sign is locally slaved to the 
d-component and varies randomly in space. 

We would like to note that both d-wave and s-wave superconducting phase are expected 
to exhibit glassy behavior associated with rare regions where the uniform phased order is 
strongly frustrated. We consider the existence of a glass phase which spontaneously breaks 
time reversal symmetry likely but not proven, so we have not shown it in the phase diagram 
in Fig. 3. 

V. DISCUSSION 

In this article we have focussed attention on several systems in which, at the point 
of the quantum superconductor-metal transition, the conductivity is still large compared 
to the quantum of conductance, and hence localization effects are unimportant. The key 
general aspects of this transition are: 1) The T = transition occurs at a point at which the 
superconducting order parameter is small in most of the sample, other than in a dilute set of 
locally superconducting puddles. 2) The transition is triggered by the quantum fluctuations 
of the phase of the order parameter on these puddles, whose quantum dynamics can either 
be governed by the dynamics of the Cooper instability (when the optimal puddles are small) 
or the electric field fluctuations (when the optimal puddles are large). 3) While there is 
probably a small quantum critical regime at exceedingly low temperatures and very close 
to the quantum critical point where the physics is universal and can be described by a 
suitable quantum critical scaling theory, there is a much larger regime where quantum 
phase fluctuations dominate the physics, but the long-distance properties are more properly 
described by the thermally truncated percolation of phase coherence between puddles. 

We have ignored the fundamental, but for our purposes purely academic question of 
electron localization in highly conducting samples. However, there remains the issue of the 
ultimate fate of the "metallic phase" in the true T — * limit. In 3D this is, presumably, 
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not an issue, but in 2D, where all single particle states are localized, this is an important 
point of principle. For Gid ^> 1 the temperature below which interference effects are 
relevant is exponentially small, so the issue is not of practical relevance. Still even the point 
of principle is interesting, and, in our opinion, unresolved partially because calculations 
of weak localization corrections close to the point of the metal-superconductor transition, 
where the conductivity is much larger than the Drude value, remains a challenge. 



3. Previous studies of the problem 

There have, of course, been many theoretical studies of the quantum phase transition from 
the superconducting to the non-superconducting state in the presence of quenched disorder. 
Many of these studies concerned a scenario in which there is a direct superconductor to 
insulator transition. Under circumstances in which, near criticality, kpl ~ 1, it is conceivable 
that there is a direct superconductor to insulator transition. Moreover, in common with 
the superconductor to metal transition considered by us, near such a superconductor to 
insulator transition the electron wave functions become strongly non-uniform, and exhibit 
fractal features [49I ]. Therefore it is not surprising that the order parameter at the point 
of the transition also non-uniform 
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51] . Despite the similarity of this aspect of the two 
transitions, the superconductor to insulator transition differs from the superconductor to 
metal transition in significant ways, and is outside the scope of the present paper. 

There have also been many previous studies of the superconductor to metal transition, 
starting with the mean-field studies of Abrikosov and Gorkov of the transition in a magnetic 



field. In this context, we would like to mention the paper Ref. [52j where a renormalization 
group approach was proposed for disordered 2D s-wave superconductors in the absence of 
magnetic field. The conclusion reached in this article bears some similarity to ours: The 
superconductivity is quenched under conditions such that Gid > 1- In the framework of 



the analysis in Ref. 



52 



the reason is that the diagrams responsible for suppression of T c 
by fluctuations of the phase of the order parameter are proportional to ln 3 (L^/7), while the 
weak localization corrections to the conductivity are only proportional to ln(/c^L<r). The 



theory presented in t 
experiments of Ref. 

physical origin as the power law in Eq. [7J However, there are significance differences in 



lis work was also supported by a comparison between theory and the 
53 . The corrections of order In 3 (Lt/ I) in Ref. [52] are of the similar 
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our analysis, the most important of which is that in the situations we have considered, the 
existence of a superconductor-metal transition is generically unconnected with interference 
effects governed by the parameter kpl. Therefore the quantum superconductor to insulator 
transition and 2D localization can be treated separately from each other. As a result, the 
effects we have considered are much larger than those considered in Ref. [52. 



Pioneering studies of the quantum d-wave superconductor (or more exotic supercon 
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which employed standard 



ductor) to metal transition were carried out in Refs. 
diagramatic techniques and the replica trick to implement the disorder averaging. Our re- 
sults differ significantly from those in Refs. jl?J, |64| , principally due to the fact that these 
earlier studies did not account for existence and crucial role of rare suprconducting puddles 
near the transition. 



A. Experiments on quantum superconductor-metal transitions 

There have been a vast number of experiments on the destruction of superconductivity 
in thin-film systems - too many for us to comment on here. For example there are many 
experiments (for a review see, for example, Ref. Q]) in which s-superconductivity in films is 
destroyed by disorder in the absence of a magnetic field, and the transition takes place at 
G e ff ~ 1. Since our theory is developed in the case G e // 1, the present theory has no 
direct relevance to these experiments. We will only discuss experiments in which G e ff ^> 1, 
and even here, only a very small subset of them. 

Before discussing the relation between our results and experiments we must address a 
question of terminology. An unambiguous distinction between metallic and insulating states 
can be made only in the limit T — > 0. The metallic state has finite resistance in this limit, 
while the insulating state has infinite resistance. (The superconducting state, of course, has 
zero resistance.) The complication concerns the way in which finite temperature data is 
extrapolated to the T — > limit. 

One relatively widely used criterion is to study the sign of the dimensionless quantity 

R T = d\og[p]/dlog[T] (60) 

at the lowest accessible temperatures, and to identify the insulating state with Rt < 0, 
a superconducting state with Rt > 0, and a metallic state with Rt ~ 0. Clearly, in the 
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extremes, this is a sensible criterion, since in order for the resistivity to either vanish or 
diverge in the T — > limit, Rt must have the stated sign. The problem comes when Rt in 
the accessible range of temperatures is relatively small in magnitude. There are certainly 
well documented ways for a metal to exhibit Rt < (for instance, in the Kondo effect), so 
an observation of Rt < cannot be safely taken as evidence that p will diverge as T — > 0. 
Of course, it is also common for metals to have Rt > 0, so this observation, by itself, 
cannot be taken as a sure indication of a zero temperature superconducting state. It is 
always also important to pay attention to the absolute magnitude of the resistance in such 
discussions. Metals at low temperatures typically have resistances smaller than the quantum 
of resistance, while insulators, at low enough temperatures, always have resistances large 
compared to this value. 

1. Transition in a perpendicular magnetic field in films with large G2D 

When highly conducting films of a superconducting metal, such as MoGe, are subjected 
to a perpendicular field, two asymptotic behaviors are expected, and observed as a function 
of H: For small enough H, the resistance tends |54j to arbitrarily small values at low T, 
which can be interpreted as a superconducting state. At large enough H, the resistance 
tends to a roughly temperature and H independent value, p = pn, which can be interpreted 
as the "normal" state (Drude) value. (pr> is expected to be almost field independent so long 
as r c l <C 1 where r c (H) is the cyclotron radius.) 



It has been observed in [54|, and latter in [55J, that at small T there is a large interval 
of magnetic fields where the resistance is independent of T, and can be up to four orders of 
magnitude smaller than pp. 

Although we have not calculated the conductance, and so cannot propose direct com- 
parisons between theory and experiment, we believe that the significant enhancement of 
the conductance takes place in the interval of magnetic fields in which somewhat isolated 
puddles of the sample have a local value of H c2 > H. (It is not completely clear, to us, 
whether the observed behaviors should be interpreted as the finite T behavior of a system 
in the range of fields, < H < H c , where the system forms a gauge glass phase in the 
T —>■ limit, or whether it should be interpreted in terms of the anomalous metallic phase 
for H > H c , where even at T = there is no global phase coherence.) In any case, to explain 
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the broad range of H over which significant superconducting fluctuations occur, we must 
assume that oh > \fG2D in Eq. [32j 

In this context we would like to mention an interesting phenomenological observation 



made recently by Steiner et al [56|]. They focussed attention on the critical value of the 
magnetic field, H*, at which Rt(H) changes sign: Rt(H) > for H < H*, Rt(H) < 
for H > H*, and Rt{H) k, for H = H*. This is often identified as the point of a SIT, 
despite the fact that, on both side of the "transition," the T dependence of the resistance 
is sometimes sufficiently weak that an unbiased extrapolation to T = would yield a finite 
result. Steiner et al found that the behaviors could be sorted into two classes. In some films, 
p(H = H*) fa h/4e 2 , and in these, a scaling collapse of the data suggestive of universal 
quantum critical phenomena can be achieved, and not only is Rt{H) < for H > H*, 
but the resistance actually grows large enough at low T that it is suggestive of a truly 
insulating phase. In other films, p(H = H*) <C h/4e 2 , and in these, resistance appears 
to approach a finite "metallic" value for H on both sides of H*, and correspondingly any 
attempt to scale the data breaks down at low T. Moreover, in these low resistance films, 
the "critical" resistance, p(H = H*), is manifestly non- universal. This analysis suggests 
that there are two possible limiting behaviors - the one in low resistance films, to which the 
present analysis is applicable, and that in the higher resistance films, which may, to some 
level of approximations, be exhibiting a superconductor to insulator transition. 



However, even in the highly conducting films of Ref. 54J, experimental indications of 
glassy behavior has not been reported for superconducting films in a perpendicular magnetic 
field, contrary to our expectations. 

2. Transition in a parallel magnetic field in films with large G2D 

Most studies of superconducting films involve relatively heavy elements, such as Mo or 
even Pb, so that the spin-orbit scattering rate is substantial and A r so > 1. However, in 
Ref. |57|, Wu and Adams studied aluminum films where A t so ^> 1. In these experiments, it 
has been observed that in the vicinity of H\\ = , the time dependence of the resistance 
exhibits long time relaxations with characteristic times of order 10 3 sec. During this period of 
time, the resistance changes by orders of magnitude and exhibits avalanche- like jumps. This 
is the characteristic dynamics of a glassy system. We think that this behavior is compatible 
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with our theory, as discussed in Section 3B. We do not know of any experiments reported 
to date on the quantum transition between this superconducting glass state and the normal 
metal - we believe such experiments could critically test the ideas presented here. 



3. Transition in d-wave superconductors as a function of disorder 

The cuprate high temperature superconductors are the best established example of a 
d-wave superconductor. Here, the critical temperature, T c , is known to vary strongly as a 
function of the doped hole concentration, x, producing two quantum critical points at which 
T c vanishes: a lower critical doping concentration, x%, on the "underdoped" side, and an 
upper critical concentration, X2, on the "overdoped" side of the phase diagram. On the 
underdoped side of the superconducting dome, the thermally accessible normal state above 
T c is manifestly not a good Fermi liquid. Moreover, with increasing underdoping, these ma- 
terials frequently appear to undergo a superconductor to insulator transition with a critical 



resistance that is typically large compared to h/4e 2 58|, |59| Thus, the present considera- 
tions may not be applicable. (However, in some instances of very high quality crystals of 
YBCO, the normal state revealed upon quenching superconductivity by underdoping can 
be somewhat metallic [6 01].) 



It is still unclear to what extent a weak-coupling, Fermi liquid based approach is valid, 
even in the "overdoped" regime of these strongly correlated materials. If we assume that, 
despite the uncertainties inherent in the strong correlation physics of the cuprates, some of 
the more robust of our findings apply to the cuprates as T c — > with overdoping, there are a 
number of interesting predictions we can make, none of which (to the best of our knowledge) 
have so far been observed experimentally. 

1) There should be a transition from a globally d-wave to a globally s-wave supercon- 
ducting state at a doping concentration, Xd- S , which is less than the critical doping, x 2 , at 
which T c vanishes. While even for Xd- S < x < X2, any local probe will see a d- wave-like 
gap structure, global phase sensitive measurements should record an s-wave state. (Some 
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evidence of such a transition may already be present in the experiments of Ref. 

2) For x near 22, the superconducting state should consist of dilute puddles in which the 
pairing is strong, floating in an otherwise metallic sea. (Indirect evidence of such a situation 



in LSCO has been presented in Ref. l62l .) 
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3) In the metallic state with x > £2, the conductivity at low temperature should diverge 
as s — ► X2, the Hall resistance should vanish, and the Wiedemann- Franz law should be 
increasingly strongly violated, in the sense that the conductivity should be greater than 
anticipated on the basis of the thermal conductivity. 

Finally we would like to mention that there are various other candidates for experimental 
studies of the above discussed effects. There is a growing consensus that there are multiple 
other examples of d-wave superconductors, including in the "115" family of heavy fermion 
superconductors and some organic superconductors. Moreover, Sr 2 Ru0 4 is known to be a 
p-wave superconductor. It has already been demonstrated [63| that superconductivity in 
these materials is suppressed by disorder when the parameter kpl of the system is still much 
larger than unity. Though the present theory was carried out specifically for the d-wave 
(i.e. spin singlet) case, we think that it is qualitatively applicable to the p-wave (i.e. spin 
triplet) case as well, at least in the presence of spin-orbit coupling. 

We acknowledge useful discussions with M. Feigelman, D. Fisher, A. Kapitulnik, D. 
Khmelnitskii, B. Shklovskii, M. Skvortsov, O. Vafek, and A. P. Young. The research was 
supported by NSF Grant No. DMR-0704151, and DOE Grant DE-AC02-76SF00515. 

VI. APPENDIX: DERIVATION OF THE EFFECTIVE ACTION 

We now sketch representative calculations required for the derivation of the effective 
actions presented in Eqs. [lj which describes the quantum dynamics of the local supercon- 
ducting "puddles" in the regime in which the mean-field solution is highly inhomogeneous. 
To begin with, we consider the effective Euclidean action, S[A], to be a functional of the pair- 
field, A, obtained by performing a Hubbard- Stratonovich transformation on an underlying 
microscopic Hamiltonian, and then integrating out the electronic degrees of freedom. Ap- 
proaching the transition from the non-superconducting side, we assume that the magnitude 
of A is everywhere small, so we can expand the action in powers of A 



where -^(r, r', t — t') is an appropriate imaginary time ordered four-fermion correlation func- 
tion (which is dependent on the precise configuration of the quenched disorder) and . . . 
represents higher order terms in powers of A. K and other response functions that enter 
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t-t')A(r',t') + ... 
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the higher order terms in the effective action are evaluated in the normal state, i.e. they 
reflect the physics of disordered metals, not the superconducting state. 
The time Fourier transform of K generically has the structure, 

K(r, r'; uo) = K (r, r') + |w|#i(r, r') + . . . (61) 

where . . . means higher order terms in powers of uo. The presence of the non-analytic 
dependence is generic in a metal, and reflects the fact that in real-time, superconducting 
fluctuations have an exponential dependence on time; they decrease exponentially if the 
normal metal state is stable, and increase exponentially if the normal metal is unstable. 

In disordered systems K(r, r'; uj) is a random function of the coordinates. Consequently 
near the point of the quantum phase transition the distribution of the order parameter can be 
visualized as a sequence of superconducting puddles separated on by a large distance. More 
precisely, at the saddle point level, a superconducting state occurs whenever K (by which 
we mean the integral operator corresponding to Kq) has at least one negative eigenvalue, 
Ko^a — £a^a- I n other words, if MinjeJ is the smallest eigenvalue of K Q , then the super- 
conducting state occurs when Min[e a ] = . Generically, the smallest eigenvalues (i.e. states 
deep in the "Lifshitz tails") are associated with wave-functions that are spatially localized 
in regions of the system that are anomalously favorable for superconductivity. However, the 
nature of these localized solutions (i.e. the spatial extent of the localized state), and the 
distribution of eigenvalues in the tails of the distribution depend on the circumstances, as 
we discuss in Sections II - IV of the paper. 

The full saddle-point value of A sp , obtained by minimizing S[A], can be expanded in 
terms of 

A sp (r) = A A(r) « £ A A(r) (62) 

a e a <0 

In this expansion, A a can be approximately interpreted as the superconducting amplitude 
on puddle a. One trouble with this, however, is that, like Wannier functions in a crystal, 
the wave-functions $ Q are not quite as localized as they should be, because they have small 
admixtures of the wave-function from neighboring puddles which are necessitated by the 
orthogonality condition, J <ir$* (r)$ a /(r) = 5 a , a '- 

We thus obtain Eq. [T] when we substitute the approximate expression in Eq. [62] into Eq. 
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|6~T1 The coefficient /?$ reflects the long-time dynamics of the order parameter, 




drdr / $*(r) J fC 1 (r,r / )$i(r / ). 



(63) 



The strength of the Josephson couplings between i-th and j-ih puddles J^- is given by the 
expression, 



By dimensional analysis (as well as explicit calculation) it is clear that a, and $ are given 
by Eq. [21 while J^- is given by Eq. [51 
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